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Abstract 

We introduce two graph polynomials and discuss their properties. 
One is a polynomial of two variables whose investigation is motivated 
by the performance analysis of the Bethe approximation of the Ising 
partition function. The other is a polynomial of one variable that 
is obtained by the specialization of the first one. It is shown that 
these polynomials satisfy deletion-contraction relations and are new 
examples of the V-function, which was introduced by Tutte (1947, 
Proc. Cambridge Philos. Soc. 43, 26-40). For these polynomials, 
we discuss the interpretations of special values and then obtain the 
bound on the number of sub-coregraphs, i.e., spanning subgraphs with 
no vertices of degree one. It is proved that the polynomial of one 
variable is equal to the monomer-dimer partition function with weights 
parameterized by that variable. The properties of the coefficients and 
the possible region of zeros are also discussed for this polynomial. 



1 Introduction and terminologies 
1.1 Introduction 

The auTi of this paper is to introduce two new graph polynomials and study 
their properties. The first one is a two- variable polynomial denoted by 
9g{/3,j) and the second one is a one-variable polynomial denoted by u;g{(3), 
which is obtained specialization of Oq- 

Partition functions studied in statistical physics have been a source of 
various graph polynomials. For example, the partition functions of the q- 
state Potts model and the bivariated random-cluster model of Fortuin and 
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Kasteleyn provide graph polynomials. They are known to be equivalent to 
the Tutte polynomial [Ij. Another example is the monomer-dimer partition 
function with uniform monomer and dimer weights, which is essentially the 
matching polynomial ^15j . 

The polynomial Oq comes from the problem of computing the Ising par- 
tition function defined as 

Z{G;J,h):= ^ exp (^'^ JeXiXj +'^hiXi^, (1) 

Xl,...,Xl^=±l e£E i£V 

e — ij 

where Jg and hi denote coupling constants and local external fields, respec- 
tively, and G = {V, E) is the underlying graph. In general, the partition 
function is computationally intractable and the Bethe approximation is a 
popular method for its approximation [2j. The approximation ratio, which 
evaluates the performance of this method, depends on the structure of the 
underlying graph. In particular, if the graph is a tree, the ratio is equal to 
one, i.e., the Bethe approximation gives the exact value of the partition func- 
tion. In principle, the approximation becomes more difficult as the nullity 
increases. In |30j . it is shown that the ratio is described by a multivari- 
ate polynomial 0g(/3,7). We derive the graph polynomial 6G{f3,j) as its 
two-variable version. 

The polynomial ujciP) is obtained from 9G{f3,j) by specializing 7 = 
2\/—l and eliminating a factor (1 — /3)I-^I~I^L We show that the poly- 
nomial coincides with the monomer-dimer partition function with weights 
parametrized by /?. In particular, for regular graphs, cj— polynomials are 
equal to the matching polynomials up to transformations. 

We discuss the properties of 6g and ug from the viewpoint of graph 
polynomials. The most important feature of these graph polynomials is the 
deletion-contraction relation: 

0g(/3, 7) = (1 - (3)0G\e{f3, 7) + /3^G/e(/3, l), 
UGW)=UJG\eW)+/3uG/eW), 

holds whenever e € E is not a loop. Note that the graph G\e is obtained 
from G by the deletion of the edge e, and the graph G/e is the result of the 
contraction of e. Furthermore, these polynomials are multiplicative: 

(^GiUG2 = 0GidG2 and UJG1UG2 = ^Gi^G2i 

where Gi U G2 is the disjoint union of Gi and G2. Graph invariants that 
satisfy the deletion-contraction relation and the multiplicative law have been 
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studied by Tutte [27j as the V-function. Our graph polynomials Oq and ojg 
are essentially examples of V-functions. 

Graph polynomials that satisfy deletion-contraction relations arise from 
a wide range of problems [UlII]. Most of them are known to be equivalent 
to the Tutte polynomial or to be obtained by its specialization, and thus, 
they have reduction formulae even for loops. Our new graph polynomials 
do not have such reduction formulae for loops and are essentially different 
from the Tutte polynomial. 

There have been few researches on specific V-functions except for those 
on the Tutte polynomial. The Tutte polynomial has attracted interest be- 
cause of its rich mathematical properties such as matroid invariance and 
connections to links [3H Hj. These properties are not shared by general 
V-functions. As described in this paper, our new V-functions also possess 
special properties, and thus, their investigation should be fruitful. 

The remainder of this paper is organized as follows. In Section 11.21 
the definitions and notations on graphs are described. Sections [21 [3l and 
m deal with the investigation of the ^-polynomial: the definition and basic 
properties of the 0-polynomial are presented in Section [21 the motivation 
for the definition is presented in Section [3l and the special values of Oq are 
discussed in Section [H Section [5] deals with the investigation of ujg including 
a study on the special value, P = 1- 

1.2 Basic terminologies and definitions 

Let G = {V, E) be a finite graph, where V is the set of vertices and E, 
the set of undirected edges. In this paper, a graph implies a multigraph, in 
which loops and multiple edges are allowed. A subset s of is identified 
with the spanning subgraph {V, s) of G unless otherwise stated. 

The notation e = ij is used to indicate that vertices i and j are the 
endpoints of e. The number of ends of edges connecting to a vertex i is 
called the degree of i and denoted by di . 

The number of connected components of G is denoted by k{G). The 
nullity and the rank of G are defined by n{G) := \E\ — \V\ + k{G) and 
r{G) := \V\ — k{G), respectively. 

For an edge e G E, the graph G\e is obtained by deleting e and G/e 
is obtained by contracting e. If e is a loop, G/e is the same as G\e. The 
disjoint union of graphs Gi and G2 is denoted by Gi U G2. The graph with 
a single vertex and n loops is called the bouquet graph and denoted by Bn- 

The core of a graph G is obtained by a process of clipping vertices of 
degree one step-by-step [24]. This graph is denoted by core(G). For example. 
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the core of a forest F is the graph of k{F) vertices without edges. A graph G 
is called a coregraph ii G = core(G). In other words, a graph is a coregraph 
if and only if the degree of each vertex is not equal to one. Note that the 
core of a graph is also called the 2-core [21j and can be generalized to the 
notion of the k-core [51 [^. 

2 Two- variable graph polynomial 9 
2.1 Definition 

First, we present the definition of a graph polynomial. For the definition, 
we define a set of polynomials {/n(3;)}^o inductively by the relations 

k{x) = 1, = 0, and = + /„-i(x). (2) 

Therefore, for instance, f2{x) = l,/3(x) = x, and so on. Note that these 
polynomials are transformations of the Chebyshev polynomials of the second 
kind: fn+2{2V^z) = iV^rUniz), where C/„(cose) = ^i*^. 

Definition 1. For a given graph G, 

0g(/3,7) ■■=^f3^'^UfMs)i^) e Z[/3,7], (3) 

sCE ieV 

where di{s) is the degree of the vertex i in s. 

In Eq. ([3]), there exists a summation over all subsets of E. Recall that an 
edge set s is identified with the spanning subgraph {V,s). Since = 0, 

the subgraph s contributes to the summation only if s does not have a vertex 
of degree one. Therefore, the summation is regarded as the summation over 
all coregraphs of the forms {V,s); we call these sub-coregraphs. In relevant 
papers, such subgraphs are called generalized loops [Si [9] or closed subgraphs 

[igiis]. 

The following facts are immediate from the definition. 
Proposition 1. 

(a) 9G,uG,{(3,l) = eG,W,j)9G,{(3n)- 

(b) ^Bj/3,7) = ELoa)/2fc(7)/3'=. 

(C) 0g(/3,7) =^core(G)(/3,7). 
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Xi X2 



Figure 1: Graph Xi and X2 



Example 1. For a tree T, 0t(/3,7) = 1. For a cycle graph C„ having 
n vertices and n edges, ^c„(/3,7) = 1 + For the complete graph K4, 
Ok^OS, 7) = 1 + 4/33 + 3/3^ + 6/3^72 + /3674. For the graph Xi, as shown in 
Figure [H 6xi{f3,j) = 1 + 3/3^ + /3^7^. For the graph X2, as also shown in 
Figured! 7) = 1 + 2/3 + /32 + /3V- 

2.2 Deletion-contraction relation and expression as Tutte's 
V-function 

2.2.1 Deletion-contraction relation 

We prove the most important property of the graph polynomial, 9, called 
the deletion-contraction relation. The following formula of fnix) plays an 
important role in the proof of this relation. 

Lemma 1. ^n,m G N, 

fn+m-2{x) = fn{x)fm{x) + f n-l{x) f m-l{x) ■ 

Proof. Easily proved by induction using Eq.([2]). □ 
Theorem 1 (Deletion-contraction relation). For a non-loop edge e G E, 

0g(/3, 7) = (1 - /3)eG\e(/3, 7) + /3^G/e(/3, 7)- 

Proof. Classify subgraph s in the sum of Eq. ([3]) depending on whether s 
includes e or not. The former subgraph s B e = ij yields —/36G\e + P^G/ei 
where Lemma [1] is used with n = di and m = dj. The latter subgraph s ^ e 
yields Ocxe- □ 

2.2.2 Relation to Tutte's V-function 

In 1947 [27J, Tutte defined a class of graph invariants called the V-function. 
The definition is as follows. 
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Definition 2. Let G be the set of isomorphism classes of finite undirected 
graphs, with loops and multiple edges allowed. Let i? be a commutative 
ring. A map V : ^ —)• ii is called a V-function if it satisfies the following two 
conditions: 

(i) V(G) = V(G\e) + V{G/e) if e G ^ is not a loop, 

(ii) V(GiUG2) = V(Gi)V(G2). 

Our graph invariant 9 is essentially an example of a V-function. In the 
definition of V- functions, the coefficients of the deletion-contraction relation 
are 1, whereas those of 9 are (1 — (3) and (3. However, if we modify 9 to 

%(/3,7) := (1-/3)-I^I+I^I/3-I^I0g(/3,7), 

we obtain a V-function 9:0^ Z[/3,7,/3-\ (1 - 

In Theorem 10 of [5], Bollobas et al. have constructed non-isomorphic 
A;-connected graphs that are not distinguished by deletion-contraction in- 
variants. The result implies that these graphs have the same 0-polynomial. 

2.2.3 Alternative expression of ^-polynomial 

By successive applications of the conditions of V-function, we can reduce the 
value at any graph to the values at bouquet graphs. Therefore, we can say 
that a V-function is completely determined by its boundary condition, i.e., 
the values at the bouquet graphs. Conversely, Tutte showed in [27] that for 
an arbitrary boundary condition, there exists a V-function that satisfies it. 
More explicitly, the V-function satisfying a boundary condition {V(-Bn)}n=o 
is given by 

scEn=0 

where z„ := X]j=o Q) ^"^^ ^n(s) is the number of connected 

components of the subgraph s with nullity n. 

Note that another expansion called the spanning forest expansion of 
V{G) is described in Section 5 of [5]. 

In the case of 9, Eq. derives the following expression. Although this 
theorem is a trivial consequence of Theorem [3] proved more directly later, 
we present a proof of Theorem [2] to clarify the relation to Eq. . 

Theorem 2. 

sC-B7i=0 
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Proof. It is sufficient to verify tliat 

sCEn=0 

By comparing tlie coefficients of x'^' in (1 — jEf)" = (1 — /3)~"(— /? + x)", we 
liave 

for every < A; < n. Using tliis equality and Proposition l.(b), we see tliat 
= E (fj (-ir+''^B,(/3,7) = (1,7)^-^(1 - 

Therefore, Eq. @ reduces to Eq. ([6]). □ 

Formulae ([3|) and ^ are both represented in the sum of the subsets of 
edges; however, the terms of a subset are different. Generally, a V- function 
does not have a representation corresponding to Eq. ([3]); this representation 
is utilized in the remainderr of the paper and makes the ^-polynomial worthy 
of investigation among all V- functions. 

2.2.4 Comparison with Tutte polynomial 

The most famous example of a V-function is the Tutte polynomial (multi- 
plied with a trivial factor). The Tutte polynomial is defined by 

TG{x,y) := - - (8) 

sCE 

It satisfies a deletion-contraction relation 

xTG\eix, y) if e is a bridge, 

Tg{x, y) = I yTQ\^{x, y) if e is a loop, 

TG\e{x, y) + Tc/eix, y) otherwise. 

It is easy to see that fcix^y) := (x - l)'=(<^)rG(x, y) is a V-function to 
Z[x, y]. For bouquet graphs, Tb„{x, y) = {x — In the case of the Tutte 
polynomial, Eq. (jH derives Eq. ([8]). 
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Moreover, the Tutte polynomial T is known to be matroidal, i.e., if Gi 
and G2 give the same cycle matroid, then = Tq^ holds [31]. Because 
Bn+m and Bn U Bm give the same cycle matroid, the relation 

TB„+m = TeJ^b^ (9) 

is a consequence of the invariance. Strictly speaking, T is not matroidal; 
however, it satisfies Eq. Q up to the easy factor and is almost matroidal. 

The V- functions 6 and T are essentially different. One intuitive under- 
standing is that Ob^-i shown in Proposition [TJ(b), do not satisfy Eq. Q even 
if multiplied with an appropriate factor. (If we set 7 = 0, this is not the 
case. See Proposition [2j) In the following remark, we formally state the 
difference irrespective of transforms between (/3,7) and {x,y). 

Remark. For any field K, inclusions (pi : Z[/3, 7, /S^"*^, (1 — /3)~^] ^ K, and 
(p2 '■ 'L[x,y] ^ K, we have 

(pio 6 ^ (f)2oT. 

Proof. It is easy to see that 02(^!b„)/<?!'2(Tbo) = (l)2{yT and (/'i(^b„)/0i(^So) = 
</>i(l-/3)-"0i(ELo (fc)/2fc(7)/3')- If 01°^ = 02of, then a„ := ELo {I)f2k{l')l3' 
= for some z & K, where 7' = 4>i{'y) and /3' = (j)i{/3). The equa- 
tion = a2 gives 7'^/3'^ = 0. This is a contradiction because /3 7^ and 
7 7^ 0. □ 

3 Motivation for definition 

In this section, we explain the motivation for considering the graph polyno- 
mial 9g, that is, the relation to the Ising partition function and its Bethe 
approximation. 

3.1 Definition of weigiited grapii version of ^-polynomial 

We consider the multi-variable version of 9g by attaching weights to the 
vertices and edges of G by 7 = {'~fi)i£V and /3 = {fie)e<^E respectively. Such 
a graph is called a weighted graph. We assume that the weights are real 
numbers. 

Definition 3. Let /3 = {/3e)eeE and 7 = {'yi)ii:v be the weights of G. 

sCEe^s i£V 
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If all vertex and edge weights are set to be the same, 0g(/3,7) reduces 
to ^g(/3,7). It is trivial by definition that 



&G^uGA(3n) = eGi(/3,7)eG.(/3,7), (10) 

eBo(/3,7) = i, (11) 

GG(/3,7) = 0core(G)(/3,7)- (12) 

In this definition, Qq is represented in the form of the edge states sum; 
however, it is also possible to represent it in the following form of a vertex 
state sum. This formula is important to show the relation to the Bethe 
approximation of the Ising partition function because the partition function 
is also given in the form of a vertex state sum. 

Lemma 2. 

eG(/3, - ^7%ev) = n (1 + ^^^i/3ee-"'C7"0 n 7TFT- 

xi,...,xjv=±l eeE i(zv + 

e — ij 

(13) 

Proof. From Eq. ([2]), we can easily verify by induction that 



/n(e-r') 



If we expand the product with respect to E in the right-hand side of Eq. (jl3p . 
it is equal to 



sCE ees ieV Xi=±l 

Then, the assertion follows immediately. □ 



3.2 Relation to Bethe approximation 

We demonstrate that the value describes the discrepancy between the 
true partition function of the Ising model and its Bethe approximation. A 
more detailed discussion of the same is found in [30] . 

The Bethe approximation is a method used for approximating partition 
functions of various statistical mechanical models p]. Here, we state it 
for the case of the Ising partition function. Recall that the Ising partition 
function on G for given J = {Je)e£E and h = {hi)i^v is defined by Eq. ([TJ. 
We write ilJij{xi,Xj) = exp{JijXiXj) and ipi{xi) = exp(/ijXi). 
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Definition 4. A set of functions {be{xi,Xj)}e£E and {bi{xi)}ii^v is called a 
belief [33J if it satisfies 

''^^be{xi, Xj) = bi{xi) for all i & V, Xi £ {±1}, and e = ij £ E, (14) 
be{xi,Xj) = 1 for all e = ij £ E, (15) 



n n ^^(^^) n Mxi,^^) n v^^c^;*)- (le) 

eeS ■''^ ■''^ ieV e£E i£V 

Then, the Bethe approximation of the partition function Zb is defined by the 
proportionality constant of Eq. ([E]): Zb HeeE 5757 Hi&v = IleeE V'e Hiev V'i- 

For given J and h, we can obtain a belief by an algorithm called belief 
propagation \20\ I33j. In practical situations, the algorithm stops in a rea- 
sonable time. Therefore, the Bethe approximation of the partition function 
is used in many applications [T7] . 

We show that Gg(/9,7) is equal to Z/Zb- We choose variables /3e and 
to parameterize {be{xi,Xj)}ei^E and {bi{xi)}i(=v, which satisfy Eqs. (fH|) 
and ([IS]): 

Prom the definition of Zb and Lemma [21 we see that 

X1,...,XJV = ±1 eSE ~'~ ''^ 

= eG(/3,7), 

where 7^ := — This equation implies that the approximation ratio is 



captured by the value of 6^. If the graph is a tree, we see from Eq. ([TTl) 
and (|12p that 0^ = 1, i.e., the Bethe approximation gives the exact value 
of the partition function. If the weights (3 and 7 are sufficiently small, we 
see that Qq ^ 1) i-e., the Bethe approximation is a good approximation. 
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The definition of Qq imphes that we can expand the approximation ratio 
by the sum of sub-coregraphs [HI |9l [3D] . This expansion often improves the 
approximation if we sum up some of the terms |14j . 

3.3 Transform of Ising partition function 

In the fohowing, we give the exphcit transform from (/9,7) to {J,h). 
We can always choose Ai, B^jh^, he,i, and Jg to satisfy 

^ ' -1 = exp(/i.a;i), 

+ 

1 + XjXj/JeCj"'^'^"'^^ = B~^ exp{JeXiXj + he^iXi + hejXj). 

Therefore, setting hi := h[ + Yl^^ he^i, we have 

ZiG; J, h) = 11^11 Be eG(/3, - C;')iev). (17) 

This fact shows that 0g(/3,7) gives the Ising partition function with (J, h), 
which is computed from (/9,7) as above. 

If = 1, or 7j = for all i G V, Eq. (jl7p reduces to the well-known 
expansion given by van der Waerden [28^ [3T] , 

Z{G; J, 0) = 2l^l Yl cosh( Je) ^ JJ tanh( Jg), (18) 

where £ is the set of Eulerian subgraphs, i.e., subgraphs in which all vertex 
degrees are even. This fact is deduced from /n(0) = 1 if n is even and 
/n(0) = if n is odd. 

It is well known by statistical physicists that Eq. ([T8]) can be extended 
to the following expression [lOj : 

Z{G;J,h) = 2l^l Yl cosh(Jg) ^ JJtanh(Jg) JJ cosh(/ii) ]^ sinh(/ii), 

eeE sCBees ieV;(s) i&Vo{s) 

(19) 

where Vg(s) (resp. Vo{s)) is the set of vertices of even (resp. odd) degree in 
s. Although both Eqs. pT|) and (fT9]) are extensions of Eq. (fTSj) and give edge 
subset expansions, they are different. An obvious difference is that only the 
sub-coregraphs contribute to the expansion in Eq. pT|) . 

Based on Eq. (fTTl) . we can say that the graph polynomial 0g'(/3,7) is a 
transformed Ising partition function with uniform coupling constants and 
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non-uniform external fields. In contrast, a bivariate graph polynomial in- 
vestigated in [1] is based on Eq. (jl9p . This polynomial corresponds to the 
Ising partition function with uniform coupling constants and external fields. 
A similar type of expression is also considered in |16j . 

3.4 Additional remarks on weighted graph version 

In this subsection, we present additional remarks on Qq by comparing it 
with 9q. The deletion-contraction relation given by Theorem [1] is generalized 
to weighted graphs as follows. If the weights (/3,7) on G satisfy 7j = for 
a non-loop edge e = ij, the weights on G\e and G/e are naturally induced 
and denoted by (/3',7') and {(3",f"), respectively. On G/e, the weight on 
the new vertex, which is the fusion of i and j, is set to be 7^. Under these 
conditions, we have 

QciPn) = {l-Pe)eG\eifi'n')+(3e@G/ei(3"n"), (20) 

which is proved in the same manner as Theorem [TJ 

If we set all vertex weights ji to be equal, the generalization of Theorem 
[2] holds. We write QciPi {n = l)i<^v) by 0g(/9;7) for simplicity. 

Theorem 3. 

eG(/3,7)=En^^"(i'^)^"^'^n/3- n (^i) 

sC-Bn=0 ees e&E\s 

Proof. In this proof, the right-hand side of Eq. ()2ip is denoted by Gg'(/3, 7). 
First, we check that Qq and Gg equal at the bouquet graphs. 

0B„(/3,7)=E^%i(i'^)n/5^ n (1-/5^) 
=EE('i')«^'nAEn(-ft) 
= EEl:('i.')M7)(-i)i"'-i"nft 

uCE siZuk=0 ^ ^ eSM 

= EEE(':')Q/-h)(-i)'"'-'nft- 

uCE 1=0 k=0 ^ / \ / eeu 
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Using the equality Yl]=k (") = which is obtained at /? = 

of Eq. d?!) , we have 

uCE een 

Second, we see that Qg{P-,i) satisfies the deletion-contraction relation 

eG(/3, 7) = (1 - fie)^G\e{(i\ l) + ^eQc/eiP", l) 

for all non-loop edges e, because the subsets including e amount to l3e@Q/f,{P, 7) 
and the other subsets amount to (1 — /3e)0G\e(/35 7)- 

By applying this form of deletion-contraction relations to both 0^ and 
0G, we can reduce the values at G to those of disjoint unions of bouquet 
graphs. Therefore, we conclude that = ©G- D 

A coloured graph is a graph with a map from the edges to a set of 
colours. If it is the set of real numbers, the term "weighted" is preferred. We 
can generalize the definition of V-functions to coloured graphs by allowing 
the coefficients of the deletion-contraction relation to depend on colours. 
Since 0g'(/9,7) satisfies Eqs. ([TO]) and ([20]) . it is a V-function of (edge) 
weighted graphs. An expansion similar to Theorem [3] holds for any coloured 
V-function because the proof only uses Eqs. (fTOj) and (pO|) . 

With regard to the Tutte polynomial, numerous works have focused on 
its extensions to edge- weighted or coloured versions. In [6j, the "universal" 
Tutte polynomial is constructed on coloured graphs by generalizing the or- 
dinary Tutte polynomial to the greatest extent possible. The "universal" 
Tutte polynomial derives other extensions of the Tutte polynomial such as 
the dichromatic polynomial for edge-weighted graphs given by Traldi [26] 
and the random- cluster model by given by Fortuin and Kasteleyn [I2j . 

Our extension, 0g'(/3,7), for weighted graphs resembles the random- 
cluster model defined by 

i?G(/3,K)= j^K'^^n/^- n (1-/5'=) 

sCE ees e&E\s 

because of Eq. (I2ip . The random-cluster model satisfies a deletion-contraction 
relation of the form 

Rg{P, = (1 - /3e)i?G\e(/3', + /3ei?G/e(/9", for all 6 G E. 

Note that this relation holds for loops in contrast to Gg(/9,7) as Rg{P-,k) 
is an extension of the Tutte polynomial. This difference arises from that of 
the coefficients of subgraphs s: k^^^^ and H (1; t)*"^'^^- 
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4 Additional properties of 9 and its implications 



4.1 Special values 

4.1.1 7 = case 

As suggested in Section [2231 if we set 7 = 0, the polynomial 9g{/3,0) is 
included in the Tutte polynomial. 

Proposition 2. 

M/3,0) = (l-/3r(«y(«)Tc(i 
Proof. From Proposition [TJ(b) and /2fc(0) = 1, we have 

0B„(/3,O) = (1 - /3)i-/3-' E (t)^' = (1 - + 

We also have fsji = (/3-i - Therefore, ^i3„(/3,0) = 

"^Bni^, T^)- Because V- functions are determined by the values at the bou- 
quet graphs, 0) = Tg{^, T^) holds for any graph G. □ 

This result is natural from the viewpoint of the Ising partition function. 
The Tutte polynomial is equivalent to the partition function of the q-Potts 
model [3]; if we set g = 2, it becomes the Ising partition function (with 
uniform coupling constants J and without external fields). In terms of 
the Tutte polynomial, such points correspond to the parameters (x, y) = 

j^), and thus, ?g(^, j^) is essentially the Ising partition function of 
that type. On the other hand, as discussed in Section [231 ^g(/3)0) is also 
essentially equal to the Ising partition function of that type. Therefore they 
must be equal up to some easy factor. 

We can say that the Tutte polynomial is an extension of the Ising parti- 
tion function (with uniform coupling constants and without external fields) 
to the g-state model whereas the ^-polynomial is an extension of it to a 
model with specific forms of local external fields. 

4.1.2 (3 = 1 case 

At (3 = 1, ^g(1, 7) is determined by the nullity and the number of connected 
components of the graph. 
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Lemma 3. For a connected graph G, 

Proof. We use the right-hand side of Lemma [2l which gives an alternative 
representation of Oq- If Xi / Xj, then 1 + XiXj^~^'^~^J = 0. Thus, only two 
terms of xi = • • • = xat = 1 and xi = -- - = XAr = — 1 contribute to the sum, 
because G is connected. □ 

If ^ = li^, then ^ - = 1. From Eq. ([22]), we see that 

MM)Ji^l . (23) 



Setting = 1, we also deduce from Eq. (|22|) that 

0g(1,O) = 2"(^). (24) 

4.2 Number of sub-coregraphs 
4.2.1 Bounds 

For a given graph G, let C{G) := {s; s C E, {V, s) is a coregraph.} be the set 
of sub-coregraphs of G. In the following theorem, the values ([23]) and (p^ 
are used to bound the number of sub-coregraphs. 

Although the following upper bound is proved in [30], here, we present 
both the proofs of the bounds for completeness. 

Theorem 4. For a connected graph G, 

n(G)-l 



2n(G) < |^((^)| < 



5- \/5 




The lower hound is attained if and only if core(G) is a subdivision of a 
bouquet graph, and the upper bound is attained if and only if core(G) is a 
subdivision of a 3-regular graph or G is a tree. 

Note that a subdivision of a graph G is a graph that is obtained by 
adding vertices of degree 2 on edges. 

Proof. It is sufficient to consider the case in which G is a coregraph and 
does not have vertices of degree 2, because the operations of taking core 
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and subdivision do not essentially change the nullity and the set of sub- 
coregraphs. 

From the definition of Eq. ([3]) , we can write 

6'g(1,7) = ^w{s;-f), 

sec 

where w{s;^) = HieV /rfi(s)(7)- -^^^ s G C, we claim that 

w{s;0) <1 <w{s;l). (26) 

The left inequality of Eq. (j26p is immediate from the fact that /n(0) = 1 
if n is even and /n(0) = if n is odd. The equality holds if and only if 
all vertices have even degree in s. Because /n(l) > 1 for all n > 4 and 
/2(1) = /3(1) = 1, we have w{s; 1) > 1. The equality holds if and only if 
di{s) < 3 for all i £ V. Then, the inequalities in Eq. ([2^ are proved. The 
upper bound is attained if and only if G is a 3-regular graph or Bq. For the 
equality condition of the lower bound, it is sufficient to prove the following 
claim. 

Claim. Let G be a connected graph, and assume that the degree of every 
vertex is at least 3 and di{s) is even for every i £ V and s £ C. Then, G is 
a bouquet graph. 

If G is not a bouquet graph, there exists a non-loop edge e = iojo. Then, 
E and E\e are sub-coregraphs of G. Thus, di^{E) or dig{E\e) = dig{E) — 1 
is odd. This is a contradiction. □ 

4.2.2 Number of sub-coregraphs in 3-regular graphs 

If the core of a graph is a subdivision of a 3-regular graph, we obtain more 
information on the number of specific types of sub-coregraphs. 
We can rewrite Lemma [3] as follows. 

Lemma 4. Let G be connected and not a tree. Then, we have 

n{G)-l 

1=0 

where Cn,i := ELi+i it) ^^2^^) for 1 < I < n - I and Cn,o := 2". 
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Proof. First, we note that for k > 1, 

This is easily proved inductively using Eq. ([2]). Then, Lemma [3] gives 

n{G) 



M1,7)=^b„,o,(1,7) = E(T) 

I — 1 ^ ^ 



f2k{l) + foil) 



k=l 
n{G)-l n(G) 

l=Q k=l+l \ / \ / 

n(G)-l 

1=0 



□ 



Theorem 5. Let G be a connected graph and not a tree. If every vertex of 
the core(G) has a degree of at most 3, then 

Cn(G),i = |{s G s has 21 vertices of degree 3.}\ 

forO<l< n{G) - 1. 

Proof. For a sub-coregraph s, Yli^v fd.i{s){7) = 7^' if only if s has 2/ 

vertices of degree 3. □ 



5 One-variable graph polynomial uj 

In this section, we define the second graph polynomial u by setting 7 = 
2y/^. Using Eq. ([2]), it is easy to verify that fn(,2^/^) = - n). 

Therefore, 

eG{(3,2V^) = E(-/3)l^l - dris)). (27) 

sCE ieV 

An interesting point of this specialization is the relation to the monomer- 
dimer partition function with the specific form of monomer-dimer weights, 
as described in Section 15.21 
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5.1 Definition and basic properties 



From Eq. (|22p . 0g'(1, 2^/ — l) = unless all the nullities of connected compo- 
nents of G are less than 2. The following theorem asserts that ^g(/3, 
can be divided by (1 — /3)I^I~I^L We define log by dividing that factor. 



Theorem 6. 



^g(/3) := 



(1-/3)1^1-1^1 



In Eq. (j27p . ^g(/3, 2-\/— T) is given in the summation over all sub-coregraphs 
and each term is not necessarily divisible by (1 — /3)l'^l-l^l; however, if we use 
the representation in Theorem [2l each summand is divisible by the factor, 
as shown in the following theorem. Theorem [6] is a trivial consequence of 
Theorem [71 

Theorem 7. 

scE n=0 

where /io(/3) := (1 - /3), /ii(/3) := 2, and /i„(/3) := forn> 2. 

Proof. From (b) of Proposition [1] and /m(2-v/— 1) = (-v/— 1)™(1 — m), we have 



esAi, 2V^) = Q (-1)^1 -2k) = < 



1 if n = 

2 if n = 1 
if n > 2. 



Theorem [2] gives 

sCEn=0 
sC-Bn=0 

Then, the assertion is proved. 



□ 



Example 2. 

For a tree T, cjtC/?) = 1-/3. For the cycle graph C„, a;c„(/5) = 1 + /3". For 
the complete graph K4, ojk^ (/3) = 1 + 2/3 + 3/3^ + 8/3^ + 16/3"^. For the graphs 
shown in Figure H uJxAl^) = 1 + /3 + 4/3^ and a;x2(/3) = 1 + 3/3 + 4/32. 



We list the basic properties of uj below. 
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Proposition 3. 



(a) UJG1UG2W) = ^GiW)^^G2W)- 

(b) logW) = ^G\eW) + /3wG/e(/3) if e e E is not a loop. 

(c) = 1 + (2n - l)/3. 

(d) UJg{(3) =UJ^ore(G)il3)- 

(e) i^g{P) is a polynomial of degree |V^ore(G)l- leading coefficient is 

n,ci/ (di — 1) and the constant term is 1. 

(f) Xei G*-"^^ be the graph obtained by subdividing each edge to m edges. 
Then, 

(/?) = (1 + /3 + . . . + /3— i)l^l-l^la;G(/3'"). 

Proof. Assertions (a— e) are easy, (f) is proved by \Eg\ — \ Vg\ = \EQ(,n) \ — 
and ea(^,{f3,2^) = eGW"',2^). □ 

Proposition 4. // G does not have connected components of nullity 0, then 
the coefficients of cog {(3) fl'^e non-negative. 

Proof. We prove this assertion by induction on the number of edges. Assume 
that every connected component is not a tree. If G has only one edge, then 
G = Bi and the coefficients are non-negative. Let G have M(> 2) edges 
and assume that the assertion holds for graphs with at most M — 1 edges. 
It is sufficient to consider the case in which G is a connected coregraph 
because of Proposition [31(a) and (d). If all edges of G are loops, G = Bn 
for some n > 2 and the coefficients are non- negative. If G = Cm, the 
coefficients are also non-negative, as in the case of Example [2J Otherwise, 
we reduce ujg to graphs with nullity not less than 1 by applying the deletion- 
contraction relation and see that the coefficients of w^Ve and ojG/e are both 
non-negative. □ 



5.2 Relation to monomer-dimer partition function 

In the next theorem, we prove that the polynomial oogH^) is the monomer- 
dimer partition function with a specific form of weights. 

A matching of G is a set of edges such that no two edges occupy the same 
vertex. It is also called a dimer arrangement in statistical physics [15J. We 
use both terminologies. The number of edges in a matching D is denoted 
by |D|. If a matching D consists of k edges, then it is called a k-matching. 
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The vertices covered by the edges in D are denoted by [D]. The set of all 
matchings of G is denoted by V. 

The monomer-dimer partition function with edge weights /i = {pe)eGE 
and vertex weights A = (Aj)jgv' is defined as 



Proof. We show that Hg'(— /3,A) satisfies the deletion-contraction relation 
and the boundary condition of the form in Proposition[3l(c). For the bouquet 
graph Bn, D = </) is the only possible dimer arrangement, and thus, 



For a non-loop edge e = iojo, we show that the deletion-contraction relation 
is satisfied. A dimer arrangement D G 2? is classified into the following five 
types: (a) D includes e, (b) D does not include e and D covers both io and 
jo, (c) D covers io but not jo, (d) D covers jo but not zq, and (e) D covers 
neither io nor Jq. According to this classification, Eg{—/3, A) is a sum of the 
five terms A, B,C, D, and E. We see that 



hg(m,a) := e n n 



De»eeD jeV\[D] 



We write A) if all weights fj-e are set to the same fi. 



Theorem 8. Let Aj := 1 + (dj — l)/3; then, 



Hb„(-/3, A) = 1 + (2n - l)/3 = 




^ (-/3)|D|(l + (d,„-2)/3) n ^^ 




+/3 E n 



□ GB iSV\lU] 



:Ci+/3C2. 
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In the same manner, D = Di + f3D2- Similarly, 

Den iev\[D] 

= E (-/3)""(l + K,-2)/3)(l + Ko-2)/3) n 

Den isv\[D] 

+ /3 E (_/3)|D|(2 + (d,„+d,„-3)/3) n A* 

[D]?iO,[D]?jO >7^»0>^0 
=: ^1 + /3^2. 

We can straightforwardly check that 

EG\ei-l3,X') = B + Ci + Di+E^ 

and 

/3HG/e(-/3, A") = ^ + /3C2 + /3Z?2 + /3i?2, (28) 

where A' and A" are defined by the degrees of G\e and G/e, respectively. 
Note that C2+-D2 in Eq. ()28p corresponds to dimer arrangements in G/e that 
cover the new vertex formed by the contraction. This shows the deletion- 
contraction relation. □ 

Let pG{k) be the number of k-matchings of G. The matching polynomial 
ac is defined by 

acix) = {-lfpG{k)x\^\-^K 

The matching polynomial is essentially the monomer-dimer partition func- 
tion with uniform weights; if we set all vertex weights as A and all edge 
weights as /i, we have 



Therefore, for a. {q + l)-regular graph G, Theorem [8] implies 

UG{u^)=aG[^ + quju^^^. (29) 

In |18j , Nagle derives a sub-coregraph expansion of the monomer-dimer par- 
tition function with uniform weights, or matching polynomials, on regular 
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graphs. With a transformation of the variables, his expansion theorem is 
essentiahy equivalent to Eq. ([29]); Theorem [8] gives an extension of the ex- 
pansion to non-regular graphs. 

As an immediate consequence of Eq. (j29p . we remark on the symmetry 
of the coefficients of loq for regular graphs. 

Corollary 1. Let G be a {q + 1) — regular graph (q > 1) with N vertices and 
Wk be the k-th coefficient ofuci/^). Then, we have 

WN-k = WkQ^-"^^ forO<k<N. 

5.3 Zeros of a;G(/3) 

Physicists are interested in the complex zeros of partition functions, because 
it restricts the occurrence of phase transitions, i.e., discontinuity of physical 
quantities with respect to parameters such as temperature. In the limit of 
infinite size of graphs, the analyticity of the scaled log partition function on 
a complex domain is guaranteed if there exist no zeros in the domain and 
some additional conditions hold. (See [32^ I23j.) For the monomer-dimer 
partition function, Heilman and Lieb [15] show the following result. 

Theorem 9 ([15] Theorem 4.6). ///ie > for all e £ E and Re(A-,) > for 
all j G V , then Hg(/i, A) 7^ 0. The same statement is true i/Re(Aj) < for 
all j G V. 

Because our polynomial ujgW) is a monomer-dimer partition function, 
we obtain a bound of the region of complex zeros. 

Corollary 2. Let G be a graph and let dm and dM be the minimum and 
maximum degree in core(G), respectively, and assume that dm>'^- If P G C 
satisfies oJcif^) = 0, then 



dM — 1 dm — i 

Proof. Without loss of generality, we assume that G is a coregraph. Let 
/3 = |/3|e'^ satisfy oJcifi) = 0, where < 9 < 2tt and i is the imaginary 
unit. Because ^^(O) = 1 and the coefficients of ujg{(3) are not negative from 
Proposition m we have /3 / and 9^0. We see that 

a;G(/3) = Hg(-/3, A) = Hg(|/3|, ze-*''/2A)(ie-^^/2)-l^l, 

where \j = 1 + (dj - l)/3 and Reiie-'^/^ Xj) = (1 - (dj - 1)|/3|) sin |. Prom 
Theorem 121 the assertion follows. □ 
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In particular, if the graph is a (g + l)-regular graph, the roots he on a 
circle of radius 'i/q, which is also directly seen by Eq. ()29p by combining the 
well-known result on the roots of matching polynomials [15]: the zeros of 
matching polynomials are on the real interval (— 2y^, 2y^). 

5.4 Determinant sum formula 

Let T := {C C E;di{C) = or 2 for all i G V} be the set of unions of 
vertex-disjoint cycles. In this subsection, an element C G T is identified 
with the subgraph {Vc, C), where Vc := {i G V; di{C) / 0}. A graph G\C 
is given by deleting all the vertices in Vc and the edges of G that are incident 
with them. 

In this subsection, we aim to prove Theorem [TOl in which we represent 
ujg as a sum of determinants. This theorem is similar to the expansion of 
the matching polynomial by characteristic polynomials |13j : 

acix) = 2'^^^ det[x/ - Ag\cI (30) 
CeT 

where ^g\c is the adjacency matrix of G \ C and k{C), the number of 
connected components of C. 

Theorem 10. 

ojg{u^) = J2 2^'^^^ det ([/ - uAg + u^{Dg - I)] 

where Dg is the degree matrix defined by {DG)i,j '■= di6ij and - {q^^q denotes 
the restriction to the principal minor indexed by the vertices of G\C. 

Proof. For the proof, we use the result of Chernyak and Chertkov For 
given weights /i = {ne)eeE and A = (Ai)jgi/; a \V\ x \V\ matrix H is defined 
by 

H := diag(A) - ^ ^/^JI^Ae, 

where = Eij + Ej^i for e = ij and Eij is the matrix base. In our notation, 
their result implies 

EGifi, X)=Y1 ^'^""^ det H\g\c n 
ceT eec 

If we set X-i = I + {di — 1)m^ and \/— /Ue = u, then the assertion follows. □ 



(31) 
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For regular graphs, Eqs. ([30]) and ([3T]) are equivalent because of Eq. (j29|) . 

The matrix (I — uAq + v?{Dg — I)) is well known for its appearance 
in the Ihara formula of the graph zeta function [25j. The result in [29] 
shows that the Bethe free energy and the graph zeta function are intimately 
related although mathematical relations between the result and Theorem [TU] 
are unknown. 

5.5 Values at /3 = 1 

The value of uJci^) is interpreted as the number of a set constructed from G. 
For the following theorem, recall that G*-^-* is obtained by adding a vertex on 
each edge in G = {V,E). The vertices of G(2) := (y(2),£;(2)) 

are classified 

into Vq and Va, where Vq is the set of original vertices and Va is that of 
newly added ones. The set of matchings on G^^^ is denoted by T>q(2) ■ 

Theorem 11. 

u;G(l) = |{DePG(2);[D]Dyo}l- 
Proof. Prom Theorem (Tj we have 

sC-B,s = GiU'-'UGj.(^j 
n{Gj) = l for j = l...fe(s) 

where Gj is a connected component of {V,s). We construct a map F from 
{D G Pfjca); [D] D Fq} to s C ^ as 

F(D) := {e G -E; the half of e is covered by an edge in D}. 

Then, the nullity of each connected component of -F(D) is 1 and = 

Example 3. For the graph X3 in Figure O ux^i^) = ^c^i'^) = 2. The 
corresponding arrangements are also shown in Figure [2j 

Finally, we remark on the relations between the results on ujci^) obtained 
in this paper. From Proposition [3l ujg{1) satisfies 

u;g(1) = a;G'\e(l) + '^G/e(l) if e S is not a loop. 

This relation can be directly observed from the interpretation of Theorem 
\TT\ Theorem [8] gives 

a^a(l)=E(-l)"" 11'^- 
Dec i&v\[r>] 
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I I 



A AA 



Figure 2: Graph and possible arrangements on 



which can be proved from Theorem [TT] with the inclusion-exclusion principle. 
Theorem [10] gives 



We can directly prove this formula from Theorem 1111 using a type of matrix- 
tree theorem. 
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